Effect of Coulomb blockade on STM current through a granular film 
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The electron transport through an array of tunnel junctions consisting of an STM tip and a granular 
film is studied both theoretically and experimentally. When the tunnel resistance between the tip 
and a granule on the surface is much larger than those between granules, a bottleneck of the tunneling 
current is created in the array. It is shown that the period of the Coulomb staircase(CS) is given by 
the capacitance at the bottleneck. Our STM experiments on Co-Al-0 granular films show the CS 
with a single period at room temperature. This provides a new possibility for single-electron-spin- 
electronic devices at room temperature. 



Charging effects on single electron tunneling such as 
Coulomb blockade and Coulomb oscillation have at- 
tracted much interest Recent advances in nano- 
technology enable us to fabricate small tunnel junctions 
where charging effects play an essential role. The TV 
characteristics for double tunnel junction systems have 
been extensively studied, where the step-like structure 
called Coulomb staircase(CS) is observed when the resis- 
tances between the junctions are not equal 0] . For these 
asymmetric double junction systems, the junction with a 
large resistance behaves like a bottleneck of the tunnel- 
ing current and the central island is charged through the 
other junction up to the maximum charge. The tunnel- 
ing current jumps when the maximum charge changes. 
Recently, the CS has been observed by using a tip of 
scanning tunneling microscope(STM) of nanometer-size 
in highly resistive granular films as well as metal- 
droplet systems [gU^]. For a granular film, which is con- 
sidered to be an array of tunnel junctions, the observed 
CS implies that a bottleneck exists in the conducting 
paths. However, the physics behind the CS in a gran- 
ular film is not clear because it contains many granules 
with different size and the conducting path may form a 
three-dimensional network in a thick granular film. 

Bar-Sadeh et al. have studied the STM current 
through a nonmagnetic granular film, AU-AI2O3, by us- 
ing the cryogenic STM |||J^. They observed the CS at 
temperatures T = 4.2 and 78 K, and analyzed the exper- 
imental data by using a triple barrier model. They as- 
sumed that the rate for tunneling between two granules is 
small and the number of excess electrons in each granule 
is treated independently. Because of these assumptions, 
the CS was given by the superposition of two different 
periods in their model: one was determined by the tun- 
nel process between the STM tip and a granule, the other 
between another granule and the base electrode. On the 
contrary, as we will show later, the CS has a single period 
which is determined by the capacitance at the bottleneck. 

In this Letter, we study the electron transport through 
an array of tunnel junctions consisting of an STM tip and 
a granular film both theoretically and experimentally. In 
this system, we can vary the tunnel resistance between 
the tip and a granule on the surface by changing the 
distance between them. When the tunnel resistance be- 



tween the tip and a granule on the surface is much larger 
than those between granules, a bottleneck of the tunnel- 
ing current is created in the array. Theoretically, we find 
that the period of the CS is given by the capacitance of 
the bottleneck even in a thick film with many granules be- 
tween the tip and the base electrode. We present results 
from STM experiments on 10 nm- and l/zm-thick Co-Al- 
O granular films which have the CS with a single period 
at room temperature. We propose that tunnel magne- 
toresistance (TMR) oscillates with the same period as 
the CS for magnetic granular films. Our theoretical and 
experimental studies provide a new direction for single- 
electron-spin-electronic devices at room temperature. 
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FIG. 1. The system with an STM tip and a thin(thick) 
granular film is schematically shown in panel (a) ((b)). 
Filled(hollow) circles represent granules in(out of) the con- 
ducting path. The insulating matrix is indicated by the shad- 
ing. The corresponding theoretical models for (a) and (b) are 
shown in panels (c) and (d), respectively. For the thick film, 
a Bethe-lattice network of granules, where each granule has 
three neighbors, is connected to the one-dimensional array. 

Our setup is schematically shown in Figs. |l|(a) and 
|l|(b). The current flows from the STM tip to the base 
electrode through a granular film. The system with a 
thin granular film in panel (a) is modeled by the one- 
dimensional array of tunnel junctions as in panel (c). We 
will show that our experimental results for the 10 nm- 
thick film are well explained by this model with = 3. 
On the other hand, such a one-dimensional array is not 
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appropriate for a thick granular film, because, as illus- 
trated in Fig. [l](b) , the conducting paths spread and form 
a three-dimensional network as the distance from the tip 
increases. Wc model this system by a one-dimensional 
array of iV-junctions connected to a Bethe-lattice net- 
work with 3 nearest neighbors as shown in Fig. |l| (d). 
Each junction is characterized by a tunnel resistance Rj, 
capacitance Cj , and carrying charge Qj . The number of 
excess electrons in the fc-th granule is represented by n^. 

The free energy for the state characterized by the set 
of charges {rii} = {ni,n2, . . .) is given by 



2Q 



(1) 



where Qi represents the charge at the surface, ^ is the 
number of electrons supplied by the voltage source and 
i goes from 1 to N. When an electron tunnels through 
the fc-th junction, the charge Qi deviates from its initial 
value by SQf. Let us consider the energy change due to 
the single electron tunneling, E^, where the superscript 
-|-(— ) denotes the process that an electron tunnels up- 
wardfdownward) through the fc-th junction in Figs. |l|(c) 
and 0(d). From Eq. (|l|), we obtain 
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voltage y^, above which the initial state (0, . . . , 0) is un- 
stable for electron tunneling through fc-th junction, i.e., 
E+{Q,...,Q) < 0, is given by 
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The last expression in Eq. (|^) is easily obtained by con- 
sidering the equivalent network shown in the inset of Fig. 
|(a), where 1/Cl = E.<fe V^z and 1/Cr = E,>fc V^. 
The threshold voltage is given by Vr = minfcV^'. As 
pointed out by Melsen et ai, the Coulomb blockade re- 
gion increases as the number of junctions, N, increases 
]16t as shown in Figs. ||(a) and ^(b). For large TV, the 
threshold voltage is expressed in terms of the average 
capacitance (C) as Vr — {e/2{C)){N — 1) and is propor- 
tional to the film thickness, because the total capacitance 
is inversely proportional to the layer thickness. 

For a thick film with many granules between the tip 
and the base electrode, the conducting paths form a 
three-dimensional network as shown in Figs. |^(b) and 
|l|(d). The decrease of the total capacitance for a junction 
network with increasing film thickness is much weaker 
than that for a one-dimensional array. Therefore, the 
threshold voltage Vr of a thick film remains of the same 
order of the magnitude as that for a thin film. Later we 
will show that the experimental results for 1 /im-thick 
film are well explained by considering the junction net- 
work as the Bethe-lattice as shown in Fig. |^(d). 



where represents the summation along the conduct- 
ing path, and 5i_k is Kronecker's delta function. The 
deviation JQf is determined by Kirchhoff's law and is 
independent of the number of excess electrons {rii}. The 
tunneling rate is obtained by using the golden rule ^ as 



r^({n.}) = 
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By solving the master equation for the probability of 
states p{{ni}) the tunneling current through the fc-th 
junction is obtained as 



= E 



(4) 



For simplicity, we neglect the effect of residual frac- 
tional charge cotunneling spin accumula- 
tion ||l^-[r^ and level quantization [[15[ in the granules. 

We first look at the I-V characteristics for a one- 
dimensional array of tunnel junctions with a bottleneck 
of the tunneling current between the tip and a granule 
on the surface. In Figs. |^(a) and ^(b), we present the 
numerical results for tunnel junctions with TV = 2 ~ 5 
at T = 0, iV being the number of junctions. Due to 
the charging energy in each junction, single electron tun- 
neling is blocked as long as the bias voltage V is lower 
than the threshold value Vt [p6|. From Eq. (0) the bias 
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FIG. 2. The tunneling current normalized with 

Iq = where Ec-^ = e^/C'2, is plotted against the 

bias voltage at T = 0. The number of junctions are A*' = 2 ~ 5 
from left to right. The tunnel resistance ratio is taken to be 
Ri : R2 ■■ Ra ■■ R4 ■■ Rs = 1000 : 1 : 1 : 1 : 1. The capacitance 
ratio is (a) Ci : C2 : C3 : C4 : C5 = \/2 : 1 : V3 : V5 : 
(b) Ci : C2 : C3 : C4 : C5 = 1/^2 : 1 : %/3 : ^5 : ^6. The 
tunneling current with ni = as a dominant state for N — 2 
is shown in the inset of panel (b) 

The CS is classified into two types as shown in Figs, 
^(a) and ^(b). The criterion is whether the capacitance 
Ci of the bottleneck is the smallest of all of the junctions, 
i.e., whether Vr — or Vr ^ min^,^! V^. A typical CS 
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for Vt = miiifc^i is given in Fig. ||(a), where electrons 
start to accumulate at the bottleneck once the bias volt- 
age exceeds Vr- Until the accumulated electrons tunnel 
out through the bottleneck, the voltage drop caused by 
them forbids electrons to tunnel through the other junc- 
tions. Therefore, the stable state is given by (rti, 0, . . . , 0) 
and the tunneling current jumps at the bias voltage where 
the number of accumulated electrons ni changes. This 
number ni is the minimum value satisfying the conditions 

^+(ni,0,...,0) <0, E+^^{ni,0,...,0)>0. (6) 

The first condition represents the accumulated electrons 
tunneling out through the first junction. The second in- 
dicates that electrons cannot tunnel through the other 
junctions, and is rewritten as (e^/Ci)ni > eV — V^^^ . 
Therefore, the tunneling current jumps at the bias volt- 
age V given by 

V = mmVS + {e^/Ci){ni~l). (7) 

From Eq. , one can easily see that the CS has a single 
period of e/Ci. The tunneling current for each plateau of 
the CS is approximately given by / ~ er^(ni, 0, . . . , 0). 

On the other hand, if Ci is the smallest, the tunnel- 
ing current does not jump at Vt as shown in Fig. ||(b). 
For Vt < V < miufc^i V^, the tunneling current is ap- 
proximately given by / ~ erj^(0, . . . , 0). Once V exceeds 
min^:^! Vrj^, electrons start to accumulate at the bottle- 
neck and the I-V curve shows a CS. The bias voltage V 
at which the tunneling current jumps is given by Eq.(|^) 

When the bottleneck is placed at another junction 
j(l < j < N), the stable state is (0, . . . , 0, n^-i = 
—rij, rij, 0, . . . , 0) and the period of the CS is determined 
by the capacitance Cj at the bottleneck. The CS is clas- 
sified into two types in the same way as those shown in 
Figs. 11(a) and §b). However, the criterion is whether 
the capacitance Cj is the smallest or not. 

The period of CS in a thick film shown in Figs. |l](b) 
and|l](d) is also determined from Eqs. (||) and (^), as long 
as the tip is coupled to a single granule on the surface 
and the bottleneck is created between them. Therefore, 
the voltage where the tunneling current jumps is given 
by Eq. (|). The CS has a single period determined by 
the capacitance at the bottleneck, even for a thick film. 

We have performed STM experiments on Co-Al-0 
granular films. Samples with different thicknesses were 
prepared by an oxygen-reactive sputtering with a Co-Al 
alloy target; 10 nm- and 1 ^m-thick C036AI22O42 films 
consisting of Co granules embedded in an Al-oxide matrix 
were deposited on glass substrates. For the 10 nm-thick 
C036AI22O42 film, a 200 nm-thick Co-Al alloy layer was 
inserted between the Co-Al-0 granular film and the glass 
substrate as a base electrode. A conventional STM sys- 
tem was used for I-V measurements under high vacuum. 
The I-V curves were obtained by using a platinum tip at 
room temperature, and by placing the tip on a Co gran- 
ule. The tunnel resistance between granules for Co-Al-0 



granular films is estimated to be about 10^ ^ IQ^fl from 
the average diameter of granules ('^ 3 nm), the average 
intergranular distance 1 nm), and electrical resistiv- 
ity ||5[,[l0|]. On the other hand, the tunnel resistance be- 
tween the tip and a granule on the surface, is about 
10** - lO^'f^. Therefore, Ri is 10^ - lO^* times larger 
than the other tunnel resistances. 

The experimental I-V curves for a 10 nm-thick 
C036AI22O42 film are plotted in Fig. ^ (a). There exist 
two or three granules between the tip and base electrode. 
Even at room temperature, the tunneling current shows 
a clear CS with a single period. The I-V curves for a 1 
/im-thick C036AI22O42 film are plotted in Fig. ^(a). We 
also find the CS with a single period. Note that for this 
thick film, 200 to 300 Co granules exist in the direction 
perpendicular to the film plane between the tip and the 
substrate. Let us now analyze our experimental results 
by using the theory presented above. 
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FIG. 3. (a)Experimental I-V curves for a 10 nm-thick 
C036AI22O42 at room temperature. A, B, C and D refer to 
different distance between the STM tip from the surface of 
the sample. The lateral position for A and B is different from 
that for C and D. (b) Corresponding theoretical curves in a 
triple tunnel junction system at T = 300 K. The tunnel re- 
sistance at the bottleneck is taken to be Ri = 600, 700, 1300, 
and 3500 Mil for the line A, B, C and D, respectively. The 
other tunnel resistances are R2 = R3 = lMf2 and the capac- 
itances are Ci = 4.48 x lO"^*^ F, C2 = 2.13 x 10"^*^ F and 
C3 = 3.62 X 10"^^ F for all curves. 

We first examine the I-V curves for the 10 nm-thick 
film. A triple tunnel junction model with a bottleneck 
between the tip and a granule on the surface is used for 
the calculation. The calculated I-V curves are shown in 
Fig. H (b). Parameter values were chosen in the ranges 
estimated from the experiments. We find that the theo- 
retical curves explain the experimental ones very well. 

For the thick film, on the other hand, the conducting 
paths are considered to form a three-dimensional network 
inside the film; this is a more complicated system com- 
pared to the thin granular films for which the CS was 
observed so far The value of the threshold volt- 

age Vt — 0.5V is not explained by the one-dimensional 
array model, because the threshold voltage Vt is propor- 
tional to the layer thickness in this model and is about 
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100 times larger than that for 10 nm-thick film; this is in 
contrast with the experimental data. 
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FIG. 4. (a) Experimental I-V curves for a 1 /^m-thick 
C036AI22O42 film at room temperature. A, B, C and D re- 
fer to different distance of the STM tip from the surface of 
the sample. The lateral tip position for A and B is differ- 
ent from that for C and D. (b)The equivalent network for 
the electron tunneling in the one-dimensional array, where 
the Bethe-lattice network is replaced by its total capacitance 
C/2. (c)The equivalent network for the electron tunneling in 
the Bethe-lattice network. 

As mentioned before, this discrepancy can be resolved 
by considering a network of the conducting paths. We 
describe the thick film as a one-dimensional array con- 
nected to a Bethe-lattice network as shown in Fig. |l|(d). 
The bottleneck is created between the tip and a granule 
on the surface. Therefore, the stable states are given by 
(ni, 0, . . . , 0) and the tunneling current jumps when the 
number of accumulated electrons rii changes. The volt- 
age where the current jumps is given by Eq. (0) and the 
CS has a single period of e/Ci even for a thick film. The 
equivalent network to obtain the threshold voltage for k- 
th junction is shown in Figs. ^ (b) and ^(c), where 
we assume, for simplicity, the same capacitance C for all 
junctions. The key point is that for electron tunneling in 
the one-dimensional array, the Bethe-lattice network is 
replaced by its total capacitance as shown in Fig. || (b). 
The total capacitance for the Bethe-lattice is C /2 and the 
bias voltage V at which the tunneling current jumps (see 
Eq. 0) is given by ^ = {e/2C){N + 1) + (e/C)(ni - 1). 
The experimental results shown in Fig. ^(a) are consis- 
tent with our model with N ^ 2. 

The TMR in tunnel junctions with ferromagnetic elec- 
trodes [ p"8|jTg| ] and magnetic granular films ||l^ is another 
attractive topic. Recently, the TMR oscillations in asym- 
metric double tunnel junctions with ferromagnetic elec- 
trodes have been studied [|o|j2l) . The condition for TMR 
oscillations is that the tunneling current shows the CS 
and the magnetic field dependence of the tunnel resis- 
tance is not the same for all tunnel junctions. The TMR 
oscillations with the same period as the CS may be ob- 
served for magnetic granular films. Our calculation shows 
that the amplitude of the TMR oscillation is about 5% 
at room temperature for the curve A in Fig. ^(b). 

In summary, the electron transport through an array of 



tunnel junctions consisting of an STM tip and a granular 
film has been studied both theoretically and experimen- 
tally. When the tunnel resistance between the tip and 
a granule on the surface is much larger than those be- 
tween granules, a bottleneck is created in the array and 
a CS with a single period is observed in the I-V curve. 
We predicted that the period of the CS is given by the 
capacitance at the bottleneck even in a thick film with 
many granules between the tip and the base electrode. 
Our STM experiments on 10 nm- and l/j,m-thick Co-Al- 
O granular films confirmed the CS with a single period 
at room temperature. TMR oscillations for magnetic 
granular films are also predicted. Our theoretical and 
experimental studies provide a new direction for single- 
electron-spin-electronic devices at room temperature. 
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